SU{b) unification for Yukawas through SUSY 

threshold effects 

Ts. EnkhbaJ] 

The Abdus Salam ICTP, 
Strada Costiera 11, 34OI4 Trieste, Italy 



Abstract 

It is well known that the supersymmetric finite threshold effects can induce 
substantial corrections to the Standard Model fermion masses. This opens an alter- 
native possibility to correct the problematic mass ratios of the lighter generations 
within the minimal SU{5) GUT. We show that with large soft yl-terms, one can 
achieve simple unification for lighter generations without additional Higgs multi- 
plet, while having sfermions lighter than 1 TeV. The presence of such large A-terms 
will distort the sfermion mass spectrum upon running from GUT scale down to the 
electroweak scale making it distinct from the universal SUSY breaking sector, espe- 
cially in the first two generations. The implications of these splittings are studied 
in K and D meson oscillations and in rare processes — > ir^vv' and — > 7r"''z^P, 
and in the latter case the effect is found to be important. 
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1 Introduction 



The supersymmetric (SUSY) version of the Standard Model (SM), while having a sta- 
bilized Higgs mass, displays a much better unification of gauge couplings than the non 
SUSY version. The Yukawa couplings of 6-quark and r-lepton unify at a reasonably good 
level as well: a slight discrepancy of ~ 20% can be remedied by various corrections. These 
successes fail to extend to the lighter two generations. In particular, the experimentally 
determined mass ratio between the down and strange quarks is an order of magnitude 
larger than the electron to muon ratio, which are predicted to be equal to each other if 
the minimal unification is assumed at the scale of the grand unification theory (GUT) for 
all generations. 

In most GUT models this shortcoming is accounted by adding either Yukawa interac- 
tions of a new Higgs multiplet or multiplets such as 45 in SU (5) [Ij or higher dimensional 
operators [2J. In SUSY theories another possibility of correcting these wrong relations 
opens up due to the threshold effects from SUSY breaking which was reported first time 
in Ref. [3], and applied for GUT in Ref. Since then, there have been many studies 
on these effects of SUSY breaking on the fermion masses. An incomplete list is given 
in Refs. [5]-[22]. In Refs. [H [10] the importance of these corrections for b-quark in the 
large tan/? limit has been pointed out for unified models based on 50(10). In particular. 
Hall et al. in Ref. [8j showed that the loop induced Qdf'H* interaction contributes a large 
effect to the down-type quark masses due to tan/3 enhancement. By now, these SUSY 
threshold corrections are integral part of phenomenological studies [23] and every popular 
code for SUSY spectrum includes them at least for the third generation. 

In parallel to these studies there have been many theoretical models which took the ad- 
vantage of these corrections to explain the observed fermion mass hierarchies. In Ref. [12] 
a model with 4th family and horizontal gauge symmetry has been considered where it was 
shown that one can achieve unification for all Yukawas without additional Higgs represen- 
tation or higher dimensional operators. Along the similar line, in Ref. [13j models were 
proposed where the masses of the lightest charged fermions are induced purely by the 
SUSY threshold corrections through flavor violating soft masses. In Ref. [H] Babu et al. 
showed that the CKM mixing can be induced in a similar manner in left-right models. 

All these studies could be classified in the following two categories: models which (i) 
explain the fermion mass pattern by soft parameters or instead (ii) use them to achieve 
Yukawa unifications for certain GUT models. In both cases one common trend that has 
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been observed in a number of studies was to move away from a minimal choice for the 
SUSY breaking parameters for achieving acceptable fermion masses and mixings [T6] 
|22j . Diaz Cruz et al. [16j has considered the effect of large flavor violating A-terms for 
down type quarks in the first two family. They start with the minimally unified Yukawa 
couplings then correct the wrong GUT relation via the large A-term effect for the down- 
type quarks. With the form they have chosen for the A-terms the correct effective Yukawa 
couplings and the Cabibbo mixing were obtained. To do so, on the other hand, they have 
concluded that the GUT relation can be corrected with sfermions heavier than 4.4 TeV 
and tan /5 ~ 2 while gluino is lighter than TeV to be consistent with the FCNC constraints 
mostly from fi —>■ ej. 

Soon the LHC will start and probe a new energy frontier around TeV. These corrections 
become much more interesting if the sfermions are lighter and hopefully reachable at the 
LHC. Then, presence of large A-terms could be related to the observed sfermion masses. 
This will be complimentary to any new FCNC signals associated with such corrections. 
If breakthroughs happen experimentally on both sides, eventually these corrections can 
be tested or ruled out by the experiments. Also very low tan /3 < 3 seems to be excluded 
by LEP II Higgs search analysis [M] except for very small window for tan/5 < 1 in the 
case of no left-right scalar mixings which gives no radiatively induced correction to the 
fermion masses. Thus it is desirable to study the issue of the minimal unification in the 
range of moderate and large tan (3 range. 

In this paper we study large A-terms for the down-type quarks, which are not propor- 
tional to the corresponding Yukawa couplings, nevertheless have flavor diagonal form, that 
corrects the wrong GUT ratios. This choice, as we demonstrate, will escape the FCNC 
constraints even for sub TeV sfermion masses. Such large A-terms split the masses of 
down type sfermions in the first two generations, making the spectrum distinct from usual 
universal SUSY parameters which could be probed at the LHC and/or ILC. We study the 
D and K meson oscillations, rare processes D tiuu and K — > ttuu. Although the rate 
for the D-meson is found to be four order of magnitude larger than its SM prediction, is 
still far from the reach of ongoing BESIII experiment. On the other hand we have found 
a large effect in the Kaon case in some of our solutions of the SUSY threshold induced 
unification. When the precision of the branching ratio improves [25j, it could be used for 
determining whether there is a trace of such new physics through the SM global fit. 

The structure of the paper is as follows. In section 2, I give a brief review of the finite 
corrections to the fermion masses, and explain based on qualitative arguments why we 



2 



need non minimal soft parameter for the minimal unification. The section 3 is devoted 
for the details of our numerical calculations and main results. The conclusions is given 
in section 4. The relevant formulae for the SUSY threshold corrections to the Yukawa 
couplings are given in the Appendix. 

2 Radiative Corrections to Fermion Masses in the 
MSSM 

In this section, we review the SUSY threshold corrections to the fermion masses and 
highlight the qualitative features of such corrections. This will enable us to see what 
choices of soft parameters would induce the correct level of threshold effects that are 
needed for changing the wrong GUT ratio. 

2.1 The SUSY threshold corrections 

Whenever a heavy field or fields decouple from a theory they induce finite shifts in the 
parameters of the theory. In the SUSY extension of the SM, these threshold corrections 
are induced as the SUSY partners of the SM fields decouple [3]. The decouphng scale 
is believed to be around or not much higher than the electroweak symmetry breaking 
scale if SUSY is to stabilize the gauge hierarchy. In particular, the Yukawa couplings 
receive finite threshold corrections from gaugino-sfermion, Higgsino-sfermion loops The 
corresponding diagrams are depicted in Figure [1] in the case of quarks. These corrections 
are especially important due to the fact that the down-type quarks and charged leptons 
obtain a new Yukawa interaction to the up-type Higgs doublet, which, upon electroweak 
symmetry breaking, results in tan/3 enhanced corrections 

Now we elaborate on the details of these corrections in the case of quarks. The full 
expressions for these corrections are given in the Appendix, which we have used in our 
numerical analysis. The total correction for the quarks are given as follows: 

{5m,) = (5mJ) + (5mf ) + (5m^) , (1) 

where the gluino-squark loop induced correction is given by 

{5m^).. ^ -^(m^R)m^/(m|,m?c,m?) , (2) 
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The definitions of tlie mass parameters in tlie above formula are given in the Appendix 
along with the details of the remaining two corrections. The loop function behaves approx- 
imately as /(m^,m^,m^) ~ l/(2m^). To see the qualitative features, here we concentrate 
on the finite correction from the gluino-squark loop which is usually the dominant one. 
For instance, the Bino-squark induced correction has 8as/(3Q;') ~ 0.03 factor compared 
to the above correction. Here we approximate the mass eigenvalues of the squaks by their 
soft-mass parameters as m? ~ ~ m?c, ignoring the mixings as well. If we assume 
the A-terms proportional to the corresponding Yukawa couplings the correction takes the 
following form 



From this formula we observe that the effect does not decouple in the limit of large SUSY 
breaking parameters. Also, concentrating on the term proportional to tan/5, it is easy 
to realize that the induced effect could be as large as the tree level term in the large 
tan/5 limit. For example, for tan/5 ~50, we see that the enhancement overcomes the loop 
suppression factor: 



If the soft masses are universal, we see that the induced changes are flavor universal. On 
the other hand, the needed corrections to the down and strange quark masses that fix 
the wrong GUT ratio are far from universal: If we want to make these corrections in the 
quark sector, we must increase the down quark mass while decreasing the strange mass. 
Therefore we seem to come to an inevitable situation that we should depart from the 
universal soft parameters. Indeed, by scanning the flavor universal soft SUSY parameter 
space, several groups have found no solution to the wrong GUT ratio [Ml [191 122] ■ In the 
next subsection, we elaborate on this issue. 

Although numerically less significant due to weaker interactions, the same tan/5 en- 
hancement occurs to the contributions from the neutralino- and chargino- squark loops. 
While they are numerically irrelevant for the lighter two generations, the chargino-stop 
loop induced correction to the bottom mass could be substantial due to large top Yukawa 
coupling. In particular, they give sizeable contributions to the (13) and (23) quark mixings 
[9]. In a certain part of parameter space, such corrections lead to ~ tan/5^ enhancement 
for the process — > [TT], which could bring it to the present upper bounds. Fur- 

thermore, if a slight discrepancy of ~ 13-24% in the 6-r unification is accounted by these 
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Figure 1: The diagram for the finite corrections to the quark Yukawa couplings. 

corrections, the branching fraction of the hght Higgs decay to bb is ahered substantially 
compared to the SM or general two-Higgs douplet models. 



2.2 Non— minimal soft A— terms for minimal unification 

Inevitable fact about SUSY extension of the SM is that it inflates the number of free 
parameters in the theory from two dozens to over a hundred upon SUSY breaking in its 
full generality. They introduce new sources for FCNCs and CP-violations which have been 
subject of intensive research for over two decades. Experimental constraints from — K^ 
mixing, flavor violating /i — 67 and many other processes suggest that SUSY breaking 
parameters are either flavor blind or aligned with the Standard Model flavor structure 
at an extremely high degree. Indeed most of the phenomenological studies, in particular 
collider analysis, concentrate on one of the universal SUSY breaking scenarios. This is 
certainly understandable considering the enormous size of the SUSY parameter space, 
which makes any attempt of generic study impractical. Secondly there is no compelling 
theoretical reason that points to a certain part of parameter space which differs from those 
universal ones. On the other hand, as we mentioned earlier, if one wishes to correct the 
wrong GUT ratios for the light fermions using the SUSY threshold corrections, one should 
most likely depart from a universal choice for soft terms. Such flavor non-universality in 
the first two generation could be phenomenologically quite distinct from the traditional 
universal scenarios and could lead to interesting rare decays. If sfermions are discovered 
at LHC, their spectrum could confirm/rule out this scenario. 

Here first we quantify what amount of threshold corrections are needed for the minimal 
SU{b) unification using approximate expressions. Then discuss the choices of the soft 
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SUSY breaking parameters which would induce such corrections. 

If one sets the strange quark mass equal to the muon mass at GUT scale, without 
taking into account the threshold effects, its low energy value would be around ~ 200 
MeV which is greater by factor of ~4 than its experimentally measured value. To get the 
correct mass value we need corrections of order ~ 150 MeV. The approximate estimate 
of the gluino-squark induced term for s-quark, for TeV soft masses can be written as 
follows: 



5ras — — —V {As cos /5 — y^yU sin /3) ^ 



±25.6Mevf^) (^^) fi^^^f 
VO.l/ V700GeV/ V m-s J 



5.0 X 



As \ ( 10 ^ 
-500GeVy Vtan/5 



+1.2 X 



^^ \ f y. 



.iTevy vi.2 X io-2yj ■ 

The plus (minus) sign belongs to a positive (negative) value for the gluino mass parameter. 
Here we made an approximation /(m?^, m?c, m?) ~ — 1/ (2m?), where m? ~ 1 TeV^. This 
crude estimate shows that the A-term, if chosen to have a large value, in spite of tan /3 
suppression, can be quite important and could even become the dominant source of the 
threshold correction. Such large values are subject to the metastability condition, which 
will be discuss shortly. To contribute constructively with the /x-term part one should 
choose the sign for Ag to be opposite to that of /i. 

If tan/3 > 30 one approaches the stability with limit of \As\ < 1.75ms TeV (See 
Eq. ([8])). Such a large value can be easily accommodated by slightly increasing one of 
the soft masses in the condition of metastability given by Eq. ([8]). To reduce the too 
large value of the s quark mass due to the unification condition, the net effect must give 
negative contribution which can be accommodated by the following choices: (i) positive 
ntg and Ag with preferably negative /i-term, (ii) negative nig and As with preferably 
positive /i-term. This choice lowers the down quark mass which has to be increased to 
fTT'd,exp- Fortunately one can see from the discussion of the s-quark case, it is much easier 
to alter the (i-quark mass via A-term due to its tiny Yukawa coupling. For example, with 
similar choice of parameter, A^ ~ 20-30 GeV is sufficient to induce the needed correction, 
which is well within the metastability limit of Eq ([8]). If we try to remedy the lighter 
two generations only using /i-term, without relying on a large A-term, we must induce 
~ —0.35 MeV change for the d-quark. Since everything is specified by the unification, 
the only freedom left is the choice of the soft masses which translates to the ratio for their 
loop functions to be Id/Is — 4.7 (/^ = I{m^^,m'^c,'m'^)). Such a large mass splitting in 
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the first two generation is unlikely to survive severe FCNC constraints for sub TeV scalar 
masses. In any case, at low and moderate values of tan /3, the //-term part cannot induce 
enough effects. This leaves us with the choice of a non-minimal A-term for either d or 
s-quark. 

Let us consider the case where we correct the d-quark mass by a large A-term while 
use fi tan (3 for correcting the s-quark mass. To have a substantial correction for the 
s-quark one must choose a larger value for the fi tan /3-term. Such a choice would, at the 
same time, reduce the 6-quark mass by a potentially large amount. This reduction can not 
be too large, otherwise it would jeopardize the already somewhat good 6-t unification. In 
the MSSM, the RG runnings of the r and b Yukawas from Mz do not give a unified value 
at GUT. Instead one obtains mr{MauT)/'mb{MGUT) — 1-13 1.24 depending on tan/3. 
Since we will impose the unification condition mr{MGUT) / Tni){MGUT) = 1 by equating 
the Yukawa coupling of the 6-quark to that of the r-lepton, its low energy value before 
the threshold correction will be higher than the experimentally determined value. Thus, 
to bring to an agreement one must choose the value of the SUSY breaking parameters 
in such a way that the net threshold correction is around —24% to —13% depending 
on the value of tan j3. Nevertheless this is a much smaller percentage change compared 
to the lighter generations: With the same mass as muon at the GUT scale one gets a 
factor of ~ 4 bigger strange quark mass, which requires 75% reduction from the radiative 
corrections. Therefore, without the effect from large A-terms, we have to choose quite 
heavier sbottom mass as can be seen from 

((5m,/m°) /, ml 

— — - ~ — ~ — I ~ 3 ^ 5. (6) 

{dmb/m°) lb mj 

Here we again come to the conclusion that the squark masses have to be very different 
from each other. To summarize, our heuristic arguments show that when the wrong GUT 
ratio of the fermion masses are corrected by the SUSY threshold effects in the down-quark 
sector, the following options are available on their soft parameters: 

(i) One of the lighter generations has a large A-term, while the sbottom soft masses 
are heavier than the remaining generations by a factor of ~ 1.5 to 2.4. 

(ii) The A-terms are proportional to their Yukawa couplings, while the soft masses 
differ from each other by large amount. In this case the FCNC constraints require them 
to be very heavy in the range of tens of TeVs. 

(iii) The s and d quarks both have a large 74-term, where soft masses can be chosen 
to be degenerate and not very heavy. 
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Certainly any of these options are viable. We find the last one interesting due to its 
potential implication on LHC phenomenology, and choose it for our numerical study. 

2.2.1 The metastability condition for A— terms 

Beside phenomenological constraints on the soft SUSY breaking parameters there are 
indirect constraints coming from the requirement of the stability of our vacuum |26j-|28j. 
In Ref. [15] Borzumati et. al. have discussed the implication of this condition for models 
of radiative fermion masses. To be absolutely stable against decaying into color/charge 
breaking vacua along the Z)-flat direction where = \fj \ = \H\, the following condition 
must be met the for Yukawa and the trilinear A-term, Yff^fH and Aff^fH: 

\Af \ < \yf\m, 

= - (mj + mj, + mjj + /x^) (7) 

where m|, m^^ and mfj are the soft masses of the fields connected by the A-term. This 
choice of absolute stability, which has no physical justification, severely restricts the size 
of A-term. When applied to the off-diagonal entries it is more constraining than FCNC 
processes [2^] . 

On the other hand if we require only metastability, namely the age of the unstable 
vacuum is longer than the age of the Universe, the constraint from numerical analysis 
gives [30] [32]: 

^ < 1.75. (8) 

m 

This metastability condition allows much more relaxed parameter space compared to the 
condition of absolute stability. With a large A-term at one's disposal it is now much 
easier to correct light fermion masses by the SUSY threshold effects. Even for large tan (3 
such a large A-term can easily overcome 1/ tan/3 suppression compared to the /i-term 
part for the down type fermions and compete with or even dominate over tan /3-enhanced 
contributions. This is an appealing scenario for any GUT model where one does not need 
to introduce additional Higgs representations. 
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3 Implications of unification through SUSY radiative 
corrections 

In this section, we present our numerical study of the SUSY radiative corrections to the 
fermion masses. In the first subsection we describe the choice of the soft parameters and 
the numerical procedures. In the last part of the section, we discuss the electroweak scale 
sfermion spectrum and their experimental implications. 

3.1 The numerical procedures and the results 

Here we describe our numerical calculations and the results. For the input values of the 
fermion masses at Msusy =1 TeV we have used the results of an recent update on the 
running fermion masses |33] (see for earlier analysis), where the running masses are 
calculated in the case of SM and MSSM at 1 TeV. Their results at 1 TeV for the SM are: 

md = 2.50 ± 1.0 MeV, m, = 47±14MeV, = 2.43 ± 0.08 GeV 

m„ = 1.10 ±0.4 MeV, = 0.532 ± 0.074 GeV, = 150.7 ± 3.4GeV. (9) 

We take the central values of these results for our analysis. Here we note that the errors 
in the estimates given by Xing, et al. in Ref. [33] are from the Particle Data Group [35] 
which are notably larger compared to lattice QCD analysis [36]. Keeping this in mind 
we seek results that are as close to the central values as possible. As for the values at 
Msusy =500 GeV we have used 

mrf = 2.63 ± 1.05 MeV, m, = 51±14MeV, m,, = 2.53 ± 0.08 GeV 

m„ = 1.17 ±0.4 MeV, = 0.553 ± 0.074 GeV, = 153.6 ± 3.4GeV. (10) 

We have used extensively SOFTSUSY [37j, a C±± based publicly available code which 
calculates the MSSM spectrum. Initially soft masses are chosen with some universal 
values. First we perform MSSM running without the threshold corrections at SUSY 
breaking scale ~ 1 TeV, to the GUT scale, [igut — 2 x 10^^ GeV, and impose minimal 
SU{b) unification: 




(11) 
(12) 

(13) 



As = = Ad, 
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= ('""'•)« = ('"«)« = ('"!») 

After the GUT conditions in Eq. ( fTTj) are imposed we run down back to the SUSY breaking 
scale. We repeat this enough until we reach stable low scale values for the Yukawa 
couplings of the down-type quarks. These values will tell us then how much corrections 
we need from the threshold effects. 

We have done the two-loop RGE running to GUT scale and back to Msusy to deter- 
mine the discrepancy between the experimentally determined values and the values de- 
rived from the unification y^. = i/i^Igut- Since the effects on the charged lepton masses are 
minor, the initial choice for the leptonic Yukawa couplings before including the threshold 
corrections would be quite close to the full effective low energy values. Thus, the numer- 
ical entries for the unified Yukawa coupling matrix Y5, defined in Eq. ( IT2|) . are chosen 
by the leptonic Yukawa coupling matrix after running them to GUT scale. Therefore 
the tree level values of the light donw-type quark Yukawa couplings will differ from their 
observed values significantly leading to the wrong GUT ratio. The objective of our nu- 
merical study is to tackle this problem by identifying the soft parameters which give the 
needed corrections. For the d and s-quarks the discrepancies are practically independent 
of tan /? and we have 

Snid ~ 1.5 MeV and Sm, ~ -156 MeV for tan /? = 5 H- 50, Msusy = 1 TeV. 

6md ~ 1.6 MeV and 5m, ~ -163 MeV for tan /5 = 5 ^ 50, Msusy = 500 GeV(16) 

As for the 6-quark mass, the needed correction mildly depends on tan/5. The result is 
summarized in Table [H Since the contribution from the A-teim is subleading for 6-quark, 



tan P 


5 


10 


15 


20 


6m,{GeV), Msusy = ITeV 


-0.690 


-0.699 


-0.687 


-0.666 


SrubiGeV), Msusy = ^OOGeV 


-0.710 


-0.720 


-0.705 


-0.693 



Table 1: The needed corrections for the mass of 6-quark. 



most of the corrections should come from the /itan/5 part. Then, the similar change of 
~ 20% are automatically induced for the other generations. For a very large choice of Ab 
at low tan/5 one can still get large effect. 
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3.1.1 Parameter choices and the induced corrections 

Now we give the details of our choice for the numerical values of the soft parameters which 
would yield the needed threshold effects to give the correct effective masses. 

The choices of the initial values of the soft terms, although highly dependent on the 
SUSY breaking mechanism and the scale at which it is mediated, we choose universal 
scalar masses at GUT scale. We choose the following parameters as the input of the 
calculation: 

(i) At GUT scale: The trilinear A-teim for the 5-plet is chosen to be simultane- 
ously diagonalize with its corresponding Yukawa coupling, Y^, but not proportional to 
the corresponding eigenvalues upon digitalization: 

(Ag).^. = ttiSij ^ aoi/5- (17) 

The sfermion soft masses are chosen to be flavor-universal 

^io = '^%^ ^3x3, 
m| = m% X /3X3, (18) 

For minimal SU (5) unification, we do not have an immediate concern to change the up 
sector Yukawa. Therefore, we keep the trilinear A-term for the 10 to be minimal at GUT 
scale A-io = a^o^io- The gaugino masses are also chosen to be universal. 

(ii) The trilinear terms A^, /i-term, the soft highs masses m|^^ and m|^^ are chosen at 
SUSY breaking scale. We work in the basis the Yukawa matrices of the charged leptons 
and down type-quarks are diagonal. For a chosen set of m"^^ and m|, first we have 
determined /i and mi/2 that give the correction given in Table [H Samples of these are 
shown in Table [2] for various choices of tan (3. After this, we have scanned over A-terms 



tan (3 


5 


10 


15 


20 




-210 


-210 


-230 


-230 


^1, = "^10 


0.314 


0.314 


0.336 


0.336 


2 2 

rrij = mi 

di 5 


0.274 


0.274 


0.294 


0.294 



Table 2: The choices for the soft parameters masses at GUT scale for various choices of 
tan (3. The units of mi/2 is GeV while that of the soft masses is TeV^ 

for the d and s-quarks, A^ and Ag, until we obtain the values of the corresponding mass 
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tan j3 


5 


10 


15 


20 


fi (GeV) 


500 


550 


580 


850 


Ad (GeV) 


3.5 


6.4 


9.2 


16.6 


As (GeV) 


-280 


-460 


-760 


-900 


Ah (GeV) 


-900 


-950 


-800 


-228 


6md (MeV) 


1.50 


1.43 


1.55 


1.69 


6ms (GeV) 


-0.170 


-0.167 


-0.158 


-0.156 


(GeV) 


-0.730 


-0.732 


-0.697 


-1.0 



Table 3: The choices for the /i-term and relevant soft trilinear A-terms at low energy and 
the induced change to the down-type quark masses. 

corrections that are close to the ones given in Eq. f|T6l) . Such values of A-terms are given 
in Table [3] with the corresponding corrections. As we can see the desired corrections are 
obtained. Here we have tried to use as low values for the /i-term as possible, such that 
the fine tuning is minimal. For this reason, A^-term has been chosen to be large, except 
in the case of tan (3 = 20 for which the A-term does not contribute significantly without 
exceeding the stability condition of Eq. (jS]). 

3.1.2 The spectrum and its implications 

Due to the large A-terms in the second generation, Ag and A^, that we found for the 
Yukawa unification, the mass degeneracy in the sfermions of the first two families is lost 
during the RG running from the GUT scale down to the electroweak scale. If sparticles 
are eventually discovered at LHC or ILC, this feature of the mass spectra makes our 
approach experimentally distinguishable from other scenarios. At the same time, the 
induced splitting cannot be too large, otherwise could exceed the experimental constraints 
from the meson oscillations and other rare processes. The sfermion masses for the above 
choices of parameters are given in Table H] and we can see a clear splitting in the masses 
of the first two generations. On the other hand, in the present case such a splitting would 
be absent in the right up-squark sector, since A-terms in the up sector are chosen to 
be proportional to the corresponding Yukawas. This would change if there were large 
corrections in the up sector as well. 

The masses in Table |4] are calculated in the basis where the down-type Yukawa cou- 
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tan (3 


5 


10 


15 


20 


tan/3 


5 


10 


15 


20 




0.505 


0.507 


0.569 


0.565 


2 

Li 


0.280 


0.375 


0.295 


0.303 


2 

V2 


0.493 


0.475 


0.481 


0.440 


2 


0.275 


0.260 


0.253 


0.243 


2 

Qa 


0.233 


0.206 


0.276 


0.395 


2 


0.202 


0.175 


0.198 


0.290 


2 


0.454 


0.456 


0.515 


0.508 




0.339 


0.348 


0.378 


0.362 


^1 


0.431 


0.394 


0.338 


0.257 




0.328 


0.318 


0.294 


0.241 


^1 


0.177 


0.114 


0.207 


0.448 




0.184 


0.147 


0.184 


0.337 



Table 4: The soft sfermion mass parameters at low energy in units of TeV^. 



tan f3 


5 


10 


15 


20 


AMd X 10^^ GeV-i 


1.44X 10-2 


0.120 


0.59 


1.42 


AMk X 10^5 GeV-i 


1.59X 10-2 


0.105 


0.706 


1.55 


AMb X 10^5 GeV-i 


1.11 


6.15 


4.53 


2.93 


AMb, X 10^^ GeV^i 


1.82 


4.56 


10.1 


45.3 



Table 5: The mass splittings in K and B meson systems due to the SUSY effects. 

pling matrix is diagonal. Therefore, upon the electroweak symmetry breaking one must 
rotate the left-handed up-squark mass matrix by the CKM matrix. Then, because of 
the induced splitting Sum'^ = — m^_^, the Cabibbo part of the rotation will induce 
non zero (12) entry of order ~ XdurriQ in the left up-squark mass matrix. The imme- 
diate consequences are the appearance of non zero mass splittings in neutral D-meson 
and Kaon systems. The Dq-Dq oscillation is induced by gluino-up squark box diagrams 
while the Kq-Kq oscillation is by that of chargino-down squark box. Observe that there 
is no gluino box diagram for down sector at the leading order, since we have chosen the 
matrices for the A-terms and Yukawa couplings for the down sector to be simultaneously 
diagonalize exactly for this reasonto avoid this correction. We have calculated these oscil- 
lation rates and found them at the safe level when checked against the latest experimental 
results shown in Table El We have included here also the results for B and Bs meson and 
are found to be negligible due to the smallness of the (13), (23) CKM mixings and the 
chargino-stop contribution at low and moderate tan/5's. As we can see the effects are 
large and , in the example of tan (3 = 20, there is already a tension with the experimental 
result of Kq-Kq. 
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(1.57±0;f 8) X 10-14 GeV [2H] 


AMk 


(3.483 ± 0.033) x lO-^^ GeV [SS] 


AMb 


(3.337 ±0.006) X lO^^^ GeV [39] 




(1.17 ±0.008) X 10~" GeV [3B] 



Table 6: The experimental values of the mass splittings of the neutral meson. 



tan P 


5 


10 


15 


20 


Br{D+ Ti+uu) X lO^i 


0.0259 


0.128 


8.39 


0.699 


Br{K+ -> vr+z/z/) x 10" 


0.141 


3.73 


37.8 


19.3 



Table 7: The branching ratios for processes — > ti^vv and tt'^uu. 



The low energy experimental searches of rare processes in the charm and strange 
sectors - as they are affected most by the non universality - could shed light on the 
spectrum we have obtained. We have examined the rates of D and K meson decays into 
pion, neutrino and antineutrino. The process — > tt^uu has a very tiny branching ratio 
in the SM for both short and long distance at the level ~ 10-^^. Although current upper 
bound is still very poor, this will improve soon once BESIII experiment starts collecting 
data soon, which would reach to the level of ~ 10^®. For our scenario the level has been 
calculated and found to be lO^^^ as shown in Table[71 Eventhough four order of magnitude 
large than the SM prediction, the rate is beyond the reach of BESIII. Therefore, further 
experimental advances are needed regarding this channel. 

The processes tt^vv and — * vr^z/z/ can be estimated reliably in the SM. 

The branching ratio of the process tt^uu has been measured by E787 and E949 

collaborations and found to be {15.7±l\^)10~^^ while for the latter there is only upper 
bound 6.7 x 10^^. If the data is improved, they provide a clean probe to possible new 
physics. Our results for the K'^ decay are shown in Table [3 Indeed in some of our fits we 
find quite a large effects from our A-term induced splitting. For example, for tan/5 = 15, 
the result is somewhat larger than the experimental result. These are calculated taking 
into account only the SUSY contributions to see the effect other than that of the SM. Since 
we do not need a complex phase in any of the A-terms, the effect cannot be significant 
for the neutral Kaon case. Thus in this case, we expect the rate to be far from the model 
independent Grossman- Nir bound [IQ]. 
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Although we do not have an exphcit model which explains the neutrino oscillation 
phenomena, most GUT models accommodate the neutrino masses through the see-saw 
mechanism, wherein one assumes three standard-model singlet right-handed neutrinos 
with masses in the range ~ 10^ 10^^ GeV. If one of them has a large Yukawa coupling 
y^^ for y^^L^u'^Hu type of interaction, it could alter r-lepton Yukawa coupling significantly 
during the running between its mass and the GUT scale [HI |12]. The /3-function of the 
r Yukawa is given by: 

= P {yT)MSSM + Y^^^^^r myc<^^< Mgut- (19) 

The additional term increases the value of yr at GUT scale, therefore, increases yb com- 
pared to yr at low energy scale. Although we have not include this possibility in our 
analysis, we would like to point out that if we get a little lower value for yi, it might be 
still compatible with the experiment in some specific models with large y,y^. The result 
for tan (3 = 20, shown in Table U\ is one possible example at hand. In this case if the 
right-handed neutrino effect causes 8-10% upward deflection on the tau Yukawa coupling 
running, the resulting low energy 6-quark mass would come out correct. 



4 Conclusions 

In this paper, we have studied the possibilities of Yukawa unification for all generations in 
SUSY SU{5) through the finite radiative SUSY threshold corrections in the presence of 
flavor non-universal soft parameters. In particular, we have concentrated on the flavor non 
universal A-terms and their effect on the Yukawa couplings of lighter generation fermions. 
It is well known that these corrections are important and they can substantially affect the 
tree level ratios between the third generation Yukawa couplings. Choosing the A-terms 
of diagonal form with large values, especially in the down quark sector, we have shown 
that the SU{5) unification with minimal Higgs content for the Yukawa interactions for 
all generation is possible for a relatively light sub TeV sfermion spectrum. This is a 
welcome scenario at the dawn of the LHC. We have examined neutral D and K meson 
mass splittings and rare processes —>■ n^uu and n^uu. While the D^-meson 

decay has been found to be much more enhanced compared to the SM, it is still far from 
the next generation experiments. On the other hand, we find the latter decay could have 
a sizable SUSY contribution and could be probed soon. 

In the present scenario, due to the large A-term for the second generation, the down 
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sector squark masses are altered from its universal value substantially making them dis- 
tinct from the most widely studied scenarios. This fact could soon be checked at the 
LHC/ILC. When the unification is assumed at the GUT scale for the soft terms, the 
slepton spectrum will display a similar distortions as well. A collider analysis of these 
type of spectra are needed to be done thoroughly. 

One of the most important questions we have not addressed in the present work is 
the problem of the proton decay, which is known to be a challenge in the minimal SU (5) 
[l3l m] . When the SUSY corrections are used in a specific GUT models, one will face the 
question of the proton decay. One could potentially evade it by the mass splitting between 
the triplet and octet which alters the unification scale with a non-renormalizable operator 
|45j . Our primary motivation was to avoid exactly these type of operators. Although this 
seems to be a setback, at least the flavor part is not necessarily affected by such operator 
if used only for rasing the unification scale. These are currently under investigation. 

If experimental breakthroughs happen in both the collider and rare decay experiments, 
the sfermion spectrum and rare Kaon decays could point in the direction of the large A- 
term scenario we have advocated in the present work. We conclude by noting the fact 
that in the absence of any theoretical prejudice for any particular SUSY breaking, all the 
phenomenologically consistent parameter space cannot be studied in full generality due 
to its huge size. Instead, the ones favored by the unification, such as the minimal Yukawa 
unification, could direct us to the parts other than the universal ones. If SUSY is the 
solution chosen by the Nature for stabilizing the gauge hierarchy we may find out soon 
which one is the correct one. 
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A Appendix 



Her we compile the formulae we have used from various sources for our numerical calcu- 
lations. These are based on the results of Refs. [IHlllTj. The superpotential of the MSSM 
is given by 

W = Y::'z,frH^Q^+Yi:^drH,Q^+YiH^H,L,+fiH^H,. 
The soft SUSY breaking part of the Lagrangian is given by: 



= K),J2/.. + K).,r5+H.,),/i5 (20) 

+ \'magg + 2"^x2XaXa + i^'m^+xlXa + (h-c.)) ■ (21) 
With the above conventions, the interactions essential to our analysis are given as: 

+ {Ci}^^PL + C^}x^Piij X^dx + di [Ci^x^Ph + C^c^x' PRj Xa^x 

+ H {>^Th',Pl + xTh^Pl) Xt + I^ {xTh'.Pr + xTKPn) X." 

+ (h.c.)). (22) 

Here, the neutralino-fermion-sfermion and chargino-fermion-sfermion couplings are given 
by: 

iVfj/) = -v/2^72tan^H^Q;(0^)^^[/5+3-r/(0^)^,^^5 (23) 
N^xL'^ = -V292 {tan^^ {Qj - Tf) (O^)^^ + Tf {O^) J f//* 

- yr(o'')ui* (24) 



a' a 



CiS = Y;^{pS,)^U:% (25) 

= -92[0^)^UZ. (26) 

cT = yi'{os,)^u::^, (2?) 

C-xt^ = n'lo'j^Ut,, (29) 



'62 ^'J' 

C?J^^ = -92{0'i)^^U:;^s + Y^' {0l)^yZ^, (30) 
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where the sfermion mixing matrices and the mass eigenstates are defined as follows: 



dmg{m}) = UfM}U^\ L = uIJl^ + uI,^J\, x = (1-6), f = u,d,l (31) 

diag{ml) = U''MlU''\ = U^^^Di, x = {1 ^ 3). (32) 

Here the charged sfermion 6x6 and sneutrino 3x3 mass matrices are: 

M} = \ ^ and _m| + M| cos 2/3/2, (33) 



ml = m)^ + r/Fjf 2{cos2 /3, sin^ [3}/2 + M| cos 2/5 {r'j - Qj sin^ Ow) , (34) 
ml = m'j^ + YfYlv'^{cos'^ (3,sin^ (3}/2- MIcos2(3{t-I -Qfsin^ew) , (35) 
mlji = A f v{cos (3,- sin f3} / V2 - fi*Yfv{smf3, cos (3}/ V2. (36) 



The matrices and (0£, O^) diagonalize the neutralino and chargino mass matrices 
respectively as follows: 

dtagim^o) = O'^M^iO^'f, diag{m^-) = 0'^MciO'i)\ (37) 
{B,Ws,HlH'J = {0''rx^ {W-,H,), = {02)^x1, (W^", ^.-) = (Og)fe(38) 



The index a' of in the neutralino contribution formula in Eq. (l23l) takes value of 3(4) 
for T| = — ^(^)- The Higgs-sfermion couplings are given as 

X'J/^ = K + Pi cos/^ + ^PU^Ul* for / = d, /, (39) 

^.f^ = -^4 sin - f^i+3f/£ for / = /, (40) 

Xl^^ = D]i^cosP~fiY:^Ul,^,U^*, (41) 

^5"^ = -/^:,sin/3 + F4sin/5-4^/,+3[/£, (42) 

X^^) = D,^^cos/5 + F,;cos/3, (43) 

X'J;^ = -^.',sin/3, (44) 

where the D and F-term induced couplings are given in terms of the following expressions: 

K = i'^!, {{Tf - Qf sm' e^) Ui^l; + Qf sm' ^h/^/,+3^£+3) , (45) 
v2cost^t^ /a « / 

2v^cos^^H^ ^'^ ^ ^ 

^4 = (( + (^Z^/)^.^^^^^^^ (47) 

Ky = 0. (48) 
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A-' = ^(-ta„9„,K)JO")^+(0")^JO-)J, (49) 
XT = ^(tan%(0«)^,(0'')^-(0")j0«)J, (50) 
= *(0g),.,(0F),,, (51) 



The neutralino-Higgs and the chargino-Higgs couphngs are given by 

92 

= S.K).,K),. (52) 
Finally, the effective Lagrangian for the Yukawa interactions arc given by 

f=d,l 

+ {Yu + SYXuiUjH'^ + {SYu)ijUiUjHl (53) 

where the corrections, SYf and 5Yj:, have contributions from the gluino-, neutralino- and 
chargino-sfermion loops. These are written as follows: 

K).- = -I'^^t^-^^'+at/J/K,-?,,-!), (54) 

ivfLif) ivrRif)* 

/jcyJvA _ '"iax ^''jay x^^^^m a I ( m'^ rr?~ rr? .\ 

^^iax ^^jbx \1{N) t( 2 2 2 \ tr-r-\ 

- ^2 ^afe ^X°^X^^(^/.'^X°'^X?j' (55) 

(5Y^) = Xl^^m--I (m'r ,m^m^-) 

V ^ Jij IQ^2 xy Xa V dy^ Xa J 

^L{u)^R{u)* 

'^iax '^jbx \ 1(C) t( 2 2 2 \ ira\ 



167r2 

-,L(d)^i?((i). 

257^2 ^^xy ""Xa " V"""^ '""'"!/' ""Xa 

Kb 'rn-- m-- 1 m, , m , m , (57) 



^iax ^jay vKu) t ( 2 2 2 



IQ^2 ab Xa Xb y d^^ Xo ' X^ 

/^^C\ ^ _ ^»a. ^ja, X^M^ 7 (^2 2 2_\ 

V ' /i7 1fi-7r2 Xa V XaJ 



167r2 

^Lil)^R{l)* 

'^iax '^jbx \l(C) T I 2 2 2 \ /roN 

16^^'^' "^x.-"^x,-^ (^i>.' ^x.-' ^x J ' (58) 

To obtain the corrections 6Yj one should replace (X^,A^) by (X^, A^) everywhere. In 

the above formulae, if one works in the basis where the down (up) type quark Yukawa 

matrix Yd {Yu) to be diagonal then the up (down) type Yukawa matrix must be chosen 

as Yf'^av^km (^/"^KL)- The loop function, I, is given by 

^ _ xy ln{x/y) + yz ln{y/z) + zx \n{z/x) 

{x - y){y - z){z - x) ' ^''^ 
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which has the following behaviors in various limits of its arguments as follows 
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